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THE VLASOV-MAXWELL-BOLTZMANN SYSTEM IN THE 

WHOLE SPACE 

ROBERT M. STRAIN 


Abstract. The Vlasov-Maxwell-Boltzmann system is a fundamental model 
to describe the dynamics of dilute charged particles, where particles interact 
via collisions and through their self-consistent electromagnetic field. We prove 
the existence of global in time classical solutions to the Cauchy problem near 
Maxwellians. 


1. The Vlasov-Maxwell-Boltzmann System 


The Vlasov-Maxwell-Boltzmann system is a very fundamental model to describe 
the dynamics of dilute charged particles (e.g. electrons and ions): 


( 1 ) 


dtF++v V.E+ + ^{e + -xbY V„E+ = Q++(E+,E+) -K Q+-{F+,F.), 
m+ \ c / 


dtF. + v-V^F- - — (e+-xb) -V^F- ^Q—{F-,F-)+Q-+{F-,F+). 
rri- \ c J 


The initial conditions are ^±(0, a;, v) = Fo^±{x, v). Here F±(t, x,v) >0 are number 
density functions for ions (-I-) and electrons (-) respectively at time t>0, position 
X = {xi,X 2 ,X 3 ) G and velocity v = {vi,V 2 ,V 3 ) G The constants e± and m± 
are the magnitude of their charges and masses, and c is the speed of light. 

The self-consistent electromagnetic field [E{t,x), B{t,x)] in (1) is coupled with 
F{t,x,v) through the celebrated Maxwell system: 

dtE — cVa; X B = —AttJ' = — 47r / v{e+F^ — e_F_}dw, 

dtB + cWx X E = 0, 


with constraints 


Va; -5 = 0, Va; • E = 47rp = 47r 


{e+E+ — e_E_}du, 


and initial conditions E{{),x) = Eo{x), B{0,x) = Bo{x). 

Let gA{v), 9 b{v) be two number density functions for two types of particles A 
and B with masses and diameters cr^ {i G {A,B}). In this article we consider 
the Boltzmann collision operator with hard-sphere interactions [3]: 

( 2 ) Q^^{gA,gB) = f \{u - v) ■ u>\{gA{v')gBiu') - gA{v)gBiu)}duduj. 

Here co G and the post-collisional velocities are 

(3) V = V -(?; — u) ■ uj\uj, u = u -\ -(u — u) ■ uj\uj. 

ruA + ruB ruA + ms 

Then the pre-collisional velocities are v,u and vice versa. 
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Recently global in time smooth solutions to the Vlasov-Maxwell-Boltzmann sys¬ 
tem near Maxwellians were constructed by Guo [11] in the spatially periodic case. 
Then convergence to Maxwellian for large times with any polynomial decay rate was 
shown for both the Vlasov-Maxwell-Boltzmann system and the relativistic Landau- 
Maxwell system [17] by Strain and Guo [18] in the periodic box. Yet the corre¬ 
sponding questions for the Cauchy problem have remained open. 

A simpler model can be formally obtained when the speed of light is sent to 
infinity. This is the Vlasov-Poisson-Boltzmann system, where B{t,x) = 0 and 
E{t,x) = \7^(f){t,x) in (1). Here the self-consistent electric potential, (j), satisfies a 
Poisson equation. There are several results for this model. 

Renormalized solutions of large amplitude were constructed by Lions [13], and 
this was generalized to the case with boundary by Mischler [16]. The long time 
behavior of weak solutions with additional regularity assumptions are studied in 
Desvillettes and Dolbeault [4]. The question of existence of renormalized solutions 
to the Vlasov-Maxwell-Boltzmann system remains a major open problem in kinetic 
theory. 

Time decay of solutions to the linearized Vlasov-Poisson-Boltzmann system near 
Maxwellian was studied by Glassey and Strauss [7, 8]. The existence of spa¬ 
tially periodic smooth solutions to the Vlasov-Poisson-Boltzmann system with near 
Maxwellian initial data was shown by Guo [10]. For some references on the Boltz¬ 
mann equation with or without forces see [2,3,6,21]. 

Next we discuss global in time classical solutions to the Cauchy problem for 
the Vlasov-Poisson-Boltzmann system. The case of near vacnum data was solved 
by Guo [9] for some soft collision kernels. Duan, Yang and Zhu [5] have recently 
resolved the near vacuum data case for the hard sphere model and some cut-off 
potentials. The near Maxwellian case was proven by Yang, Yu and Zhao [19] with 
restrictions on either the size of the mean free path or the size of the background 
charge density. After the completion of the results in this paper we learned that 
Yang and Yu [20] have removed those restrictions and further shown a time decay 
rate to Maxwellian of The method in [19,20] makes use of techniques 

from Liu and S.-H. Yu [14] and Liu, Yang and S.-H. Yu [15] as well as the study 
in conservation laws. But there seems to be serious obstacles to generalizing the 
method of [19, 20] to the full system (1) due to difficulties involved in controlling 
the electromagnetic field, in particular B(t,x). 

It is our purpose in this article to establish the existence of global in time classical 
solutions to the Gauchy Problem for the Vlasov-Maxwell-Boltzmann system near 
global Maxwellians: 




^0 ( ^ ^ ^-m+\v\^/2KTo 

e+ \ 2'kkTq J ’ 

— ( \ ^ g-m-|j)|^/2KTo 

e_ \27r/tTo } 


Accordingly, in the rest of this section we will reformulate the problem (1) as a 
perturbation of the equilibria. 
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Since the presence of the physical constants do not create essential mathemati¬ 
cal difficulties, for notational simplicity, we normalize all constants in the Vlasov- 
Maxwell-Boltzmann system to one. We normalize the Maxwellian 

^{v) = fi+{v) = AI-(t) = 

We further normalize the collision operator (2) as 

Q'^^{9A,gB) = Q{gA,gB) = / \{u-v) ■uj\{gA{v')gB{u') - gA{v)gB{u)}duduj. 

And we define the standard perturbation f±{t,x,v) to /i as 

F± = g + VJdf±- 

Let [•,•] denote a column vector so that F{t,x,v) = [F+,F_]. We then study the 
normalized vector-valued Vlasov equation for the perturbation 

f{t,x,v) = [f+{t,x,v),f-{t,x,v)], 

which now takes the form 

(4) -f T • Va; -f q{E + VX B)- V„}/ - {E ■ v}y/JIqi + Lf = |{A • v}f + T{f, /), 

with /(O, x, v) = fo{x, v), qi = [1, -1] and q = diag(l, -1). 

For any given g = [g^,g-], the linearized collision operator in (4) is 

Lg = [L+g,L_g], 

where 

(5) L±g = -2n~^/‘^Q(y/JIg±, g) - y/Jl{g± + g^}). 

We split L in the standard way [6]: Lg = v{v)g — Kg. The collision frequency is 

(6) vlv) = / \v — u\ij,{u)duduj. 

For g = [g+,g_] and h = the nonlinear collision operator is 

T{g,h) = [T+{g,h),T_{g,h)], 

with 

r±(5, h) = g~^^‘^Q{y/Jig±, y/gh±) + g,~^^‘^Q{y/Jdg±, y/Jdh^). 

Further, the coupled Maxwell system now takes the form 

(7) dtE X B =-J = - [ Vy^{f+ - f_)dv, dtB + VxXE = 0, 

JR3 

(8) V,-E = p= [ ^if+-f-)dv, V^-B = 0, 

with the same initial data E{0,x) = Eo{x), B{0,x) = Bo{x). 

It is well known that the the linearized collision operator L is non-negative. For 
fixed (t, x), the null space of L is given by (Lemma 1): 

(9) Af = span{[v^,0], [0, {1 < i < 3). 

We define P as the orthogonal projection in to the null space Af. With 

{t,x) fixed, we decompose any function g{t,x,v) = [gA{t,x,v),g-{t,x,v)] as 

g{t, X, v) = Pg(t, X, v) + (I- P)g(t, x, v). 

Then Pg is call the hydrodynamic part of g and (I — P)g the microscopic part. 
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2. Notation and Main Results 

We shall use (•,•) to denote the standard inner product in for a pair of 
functions gi{v), g 2 {v) G The corresponding norm is written \giW = 

(ffi)ffi)- We also define a weighted inner product in as 

{gi,g2)u = {v{v)gi,g2). 

And we use \-\i, for its corresponding norm. We note that the weight (6) satisfies 

-(1 + |t|) < u(u) < c(l + |u|). 
c 

We also use || • |j to denote norms in either x R^ or R^. In other words we 
use II • II to denote an norm in R^ x R^ if the function depends on {x, v). But if 
the function only depends upon a; G Ri^ then || • || will denote the norm of L^(Ri^). 
The T^(R^ X R^) inner product is written (•,•). Further, ||g||^ = {vg,g). 

Let multi-indices a and (3 be 

oc = [aoTai,a2,a^], /? = [/?i,/32)/^s]- 
We define a high order derivative as 

If each component of a is not greater than that of d’s, we write a < a] a < a 

means a < d and |a| < |d|. We also denote by Cg. 

We next define an “Instant Energy functional” for a solution to the Vlasov- 
Maxwell-Boltzmann system, [f{t,x,v),E{t,x),B{t,x)\, to be any functional £{t) 
which satisfies the following for an absolute constant C > 0: 

y{t)< Y. ll^d/WlP+ E \WE{t),d-B{t)W <C£{t). 

|a| + |/3|<N |a|<IV 

The temporal derivatives of [/o,£’o,i3o] are defined through equations (4) and (7). 
We will always use C to denote a positive absolute constant which may change from 
line to line. 

We additionally define the “Dissipation rate” for a solution as 

v{t)^\\E{t)r+ Y ii^“/wii'+ E ii5|(i-p)/wii^ 

0<|a|<A |a| + |/3|<N 

Notice that we include the electric field E(t, x) in the dissipation and further we do 
not take any velocity derivatives of the hydrodynamic part, P/. Throughout this 
article we assume fV > 4. Our main result is as follows. 

Theorem 1. Assume that [/o,ifo,i3o] satisfies the constraint (8) initially. Let 
Fo,±{x,v) = g, + Vm/o,±( a;, w) > 0. 

There exists an instant energy functional and cq > 0 such that if 

£{0) < eo, 

then there is a unique global solution [f{t,x,v),E{t,x),B{t,x)] to the Vlasov- 
Maxwell-Boltzmann system (4) and (7) with (8) satisfying 

+ [ 'L){s)ds < £^(0). 

Jo 
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And moreover F±{t, x,v) = /i + x, v) > 0. 

In particular our solutions are for N chosen large enough. The proof of 
Theorem 1 makes use of some of the techniques developed in [11,12], Yet several 
new difficulties arise which are fundamental to the full space problem for the the 
Vlasov-Maxwell-Boltzmann system and we develop new techniques accordingly. 

We expect that the techniques developed in this paper will be useful for working 
on other kinetic equations with force terms in the whole space. In particular, these 
techniques can be used to build solutions to the relativistic Landau-Maxwell system 
[17] in the whole space. And our proof of Theorem 1 supplies another construction 
of solutions to the Cauchy problem for the Vlasov-Poisson-Boltzmann system near 
Maxwellian. 

As in [11,12], a key point in our construction is to show that the linearized 
collision operator (5) is effectively coercive for solutions of small amplitude to the 
Vlasov-Maxwell-Boltzmann system: 

Theorem 2. Let lf(t,x,v),E(t,x),B(t,x)] be a classical solution to (4) and (7) 
satisfying (8). There exists Mq > 0 such that if 

(10) E {wdvmf + wm\f + WBit)\f} < Mo, 

\a\<N 

then there are constants Jq = So{Mq) > 0 and Co > 0 so that 

E {Ldy{t),d-fit))>6o\\{i-p)fit)\\i+so E 

|a|<A 0<|a|<Ar 

where G{t) = Vx ■ b(a+ + a-)dx and a±, b are defined below in (15). 

To prove Theorem 2, we split the solution / to the Vlasov-Maxwell-Boltzmann 
system (4) into it’s hydrodynamic and microscopic parts as 

(11) / = P/ + (I-P)/. 

Separating into linear and nonlinear parts, and using L{P/} = 0 (Lemma 1), we 
can express the hydrodynamic part, P/, through the microscopic part, (I — P)/, 
up to a second order term: 

(12) {dt + v VJP/ - {E ■ v}^qi = ;((I - P)/) + h{f), 
where 

(13) ;((I-P)/) = -{dt + vWx + L}{I-P)f, 

(14) hif) = -q{E + vxB)-\/xf + ^{E-v}f + T{f,f). 

Expanding Pf = [P+/, P_/] as a linear combination of the basis in (9) yields 

(15) P±/= i^a±{t,x) + Y^^b^{t,x)v^ + c{t,x)\vf'^ 

A system of “macroscopic equations” for the coefficients a{t,x) = [a+ {t,x),a-{t,x)], 
bi(t,x) and c(t,x) can be derived from an expansion of the left side of (12) in the 
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velocity variables, using (15) to obtain [11]: 


(16) 

V xC 

tj 

-b 

tj 

II 

(17) 

dfC + d^bi 

— li hi^ 

(18) 

d^bj+d^b. 

— ^ J 

(19) 

dtbi + d"a± T Ei 


(20) 

dta± 

= It + h^ 


Here = dxy The terms on the r.h.s. of the macroscopic equations (16)-(20) 
are obtained by expanding the r.h.s. of (12) with respect to the same velocity 
variables and comparing the coefficients on the two sides. Here 
lij(t, x), x) and l^{t, x) are the coefficients of this expansion of the linear term 
(13), and hc(t,x), hi{t,x), hij{t,x), hf^{t,x) and h^{t,x) are the coefficients of the 
same expansion of the second order term (14). These terms are defined precisely 
in Section 4. 

Theorem 2 is proven via energy estimates of these macroscopic equations (16)- 
(20), which are used to estimate derivatives of a, b and c. In the periodic case [11], 
the Poincare inequality was used to estimate a, b and c in terms of their derivatives 
and the conserved quantities. Further, it was observed in [12] that a, b and c 
themselves can not be estimated directly from (16)-(20). This is why we remove 
P/ from the dissipation rate in the whole space. Additionally, due to the velocity 
gradient in the equation (4), it is important to include the velocity derivatives in 
our norms. On the other hand, in our global energy estimates (in Section 5) dpFf 
behaves almost the same as P/. We therefore take velocity derivatives only on the 
microscopic part of the solution in the dissipation. Our new weaker dissipation, 
allows the proof to work. But at the same time this lack of control over the 
hydrodynamic part and it’s velocity derivatives causes many of the new difficulties. 

We use Sobolev type inequalities to bound the nonlinear terms by derivatives of 
/ in Lemma 5. This is an important step for bounding the macroscopic equations 
by our weak dissipation in the proof of Theorem 2. Lemma 5 is also used to include 
the electric field in the dissipation in Lemma 6. 

Further, in the periodic case [11] the Poincare inequality was used to estimate the 
low order temporal derivatives of bi. Instead, we notice that there is a cancellation of 
the most dangerous terms in the macroscopic equations which allows us to estimate 
the low order temporal derivatives of bi in the whole space directly in the proof of 
Theorem 2. 

Then we prove global existence via energy estimates of the Vlasov-Maxwell- 
Boltzman system in Theorem 3. New difficulties arise in these energy estimates 
for both the terms with velocity derivatives and those without because we can not 
control the hydrodynamic part of a solution by the linear operator in the whole 
space. We develop different procedures to prove the energy estimates with the 
weaker dissipation. In particular, we estimate the nonlinear part of the collision 
operator by the weaker dissipation in Lemma 7. Further, including the electric 
field in the dissipation rate (Lemma 6) enables us to control other nonlinearities 
without velocity derivatives through Sobolev inequalities. Additionally, we split the 
solution as (11) in several places because this allows us to absorb velocity derivatives 
and exploit the exponential velocity decay of the hydrodynamic part (15) in our 
estimates. 
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This article is organized as follows: In section 3 we recall some basic estimates 
for the Boltzmann collision operator. In section 4 we establish positivity of (5) for 
solutions to the Vlasov-Maxwell-Boltzmann system with small amplitude. And in 
Section 5 we prove global existence (Theorem 3). 


3. Local Solutions 


We will make use of the following basic estimates from [1,10,11]: 

Lemma 1 ([1,11]). {Lg,h) = {g,Lh), {Lg,g) > 0. And Lg = 0 if and only if 
9 = ^9. where P was defined in (15). Furthermore there is a S > 0 such that 

{Lg,g)>5\{l-V)g\l. 

From (5) it is well known that we can write L = v — K where v is the multipli¬ 
cation operator dehned in (6) and K satisfies the following. 

Lemma 2 ([10]). ///3 ^ 0, djsn{v) is uniformly bounded. For g{v) G iLl^l(M^;M^) 
and for any g > 0, there is Cjj > 0 such that 

\df3[Kg]\l<g ^ \df3'g\l + C^\g\l 

l/3'l=l/3| 


Above and below (R^;K.^) denotes the standard L^(]R^;R^) Sobolev space 
of order |/3|. Next up are estimates for the nonlinear part of the collision operator. 

Lemma 3 ([10]). Say g{v),p{v),r{v) G C)?°(R^;R^). Then 

\{df 3 T{g^p),r)\<C kU- 


/3i+/32</3 


Moreover, 

T(g^p)rdv 


T(p^g)rdv 


< Csup \v'^r\ sup 


-|l/2 

\g{x,v)\^dv ||p||. 


Using these estimates, it is by now standard to prove existence and uniqueness of 
local-in time positive classical solutions to the Vlasov-Maxwell-Boltzmann system 
(4) and (7) with (8) and to establish the continuity of the high order norms for a 
solution as in [11, Theorem 4]. Therefore, we will not repeat the argument. 


4. Positivity of the Linearized Collision Operator 

In this section, we establish the positivity of (5) for any small solution [/, E, B] 
to the Vlasov-Maxwell-Boltzmann system (4) and (7) via a series of Lemma’s. 
First we simplify the Maxwell system from (7) and (8). By (15) 

j [Vy/h,-Vy/Jj] ■ Pfdv = 0, 

J [^/M ,-'M ■ fdv = a+- a-. 

We can therefore write the Maxwell system from (7) and (8) as 

dtE - Va; X S = - J = [ [-vy^, Vyfp\ ■ (I - P)fdv, 

Jr3 

dtB -f Va; X if = 0, 

VX • E = a+ — a-, Vx ■ B = 0. 


(21) 
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It is in this form that we will estimate the Maxwell system. 

Next we discuss the derivation of the macroscopic equations (16)-(20). We ex¬ 
pand the r.h.s. of (12) in the velocity variables using the projection (15) to obtain 


-I- {dtc + d^hi}vl -f -f d^hi}viVj + {dtbi -f (9M± =F Ei}vi | ^ 

The above represents two equations. Here = dxy For fixed this is an 

expansion of l.h.s. of (12) with respect to the following basis (1 < * < J < 3): 

(22) [ui7/1,0], [0,z;*77I], [7 m,0], [0,7/1]. 

By expanding the r.h.s. of (12) with respect to the same velocity variables and 
comparing the coefficients of each basis element in (22) we obtain the macroscopic 
equations (16)-(20). Given the full basis in (22), for fixed {t,x), lc{t, x),li{t, x), 
lij{t,x), l^-{t,x) and l^{t,x), from (16)-(20), take the form 

(23) [ l{{I-P)f)-e^{v)dv. 

JR3 


Here ?((I — P)/) is defined in (13) and the {e„(u)}, which only depend on v, are 
linear combinations of the elements of (22). Similarly he, hi, hij, and h^ from 
(16)-(20) are also of the form 

(24) [ h{f)■ e„{v)dv. 


where h{f) is defined in (14). In the next two lemmas, we will estimate the right 
side of the macroscopic equations using (21), (23) and (24). 

Lemma 4. Let a = [ao, ai,a 2 , a^], jaj < fV — 1, then for any 1 < i, J < 3, 

P%ii + ir;7 + ll5%ll + ||a%±ii + ii9“eil + irf711<c ^ ll(i-p)a'^9“/ll. 

a|<l 

This Lemma is similar to [11, Lemma 7], but we prove it for completeness. 


Proof. The estimate for ff follows from (21): 

|a“f7| <c [ |777i|(i-P)a“/|du<c'|(i-p)9“/|2. 

Jr3 


Square both sides and integrate over a; G to get the estimate in Lemma 4 for J. 

For the other terms in Lemma 4, it suffices to estimate the representation (23) 
with (13). Let |q;| < N — 1. Since iy(v) < C{1 + |r!|) and K is bounded from 
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to itself. We have by (13) 

(^J 5“I((I - P)/) • = (^Ji{dt + P • V, + L}(I - P)a“/) • en{v)dv^ , 

<C [ \en{v)\dv [ \en{v)\{\il-P)dtd‘^f\^ + \v\‘^\{l-P)V,d‘^f\‘^)dv 
Js.3 JrS 

+C [ \en{v)\dv [ \eniv)\\{v - K)(l-P)d°'f\'^dv 

Jr3 Jr3 

< C{\{1 - P)dtd‘^f\l +1(1 - P)v,a“/|2 +1(1 - P)a“/|2}. 

Here we have used i9“(I — P)f = (I — P)i9“/, and the exponential decay of e„(T). 
The estimate in Lemma 4 follows by further integrating over a; G □ 

Next, we estimate coefficients of the higher order term h{f) on the right side of 
the macroscopic equations (16)-(20). 

Lemma 5. Let (10) hold for Mq > 0. Then for 1 < i j < 3, 

^ {lia“he|| + ||5“h.|| + p“/i.,|| + ||a“h±j| + ir^^ ^ p^/l|. 

\a\<N 0<|a|<AI 

Notice that the left side of this estimate contains zero’th order derivatives but 
the right side does not. We have to tighten the estimate in this way because Pf 
(with no derivatives) is not part of the dissipation rate. We use Sobolev’s inequality 
on the second order terms to remove the terms without derivatives. 

Proof. Notice that it suffices to estimate (24) with (14). For the first term of h{f) 
in (14), we integrate by parts to get 

f d'^{q{E-\-v X B) -Vyf)} • en{v)dv = f ^ {q{E v x B)jdyJ)} • en{v)dv 

3 

= ^ C“ y d,^{{d^^E + VX • er,{v)dv 

Take the square of the above and further integrate over a: G the result is 
(25) + 

This follows from Cauchy-Schwartz and the exponential decay of Vt,e„(v). 

We will estimate (25) in two steps. First assume \a — ai\ < N/2. In this case 
we use the embedding 1F^’®(]R^) C L°°(IR^) combined with the L®(]R^) Sobolev 
inequality for gradients to obtain 

sup [ \d°‘~°‘^f{x,v)\‘^dv<C [ ||a“““i/llwi.6(R3)C^w, 

a:GR3 JR3 JR3 

<c [ ||V,5“-“VfwL2(R3)d^^. 

Therefore, when jo; — ai| < N/2 and > 4 we see that 

0<|a|<2 
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The small amplitude assumption (10) completes the estimate in this case. 

Alternatively if |q; — ai| > N/2, then |q;i| < N/2 and we use the embedding 
^ to obtain 

sup {\d°^^E\ + \d°^^B\) < C ^ {\\d^d°^^E{t)\\ + ||a“a“iB(t)||}. 

a;GR3 \a\<2 

In this case 


< C'||a“-“i/ll {||9“5“iA(t)|| + ||5“a“iB(t)||}. 


Again, the small amplitude assumption (10) completes the estimate. 
For the second term of h{f) in (14) we have 


|5“{(A • v)f} ■ eniv)dv = ^ I {|(a“^F; • u)5“-“v} ' eniv)dz 


This term is therefore easily treated by the last argument. 

For the third term of /i(/), we have 

ya“r(/,/).e„(u)du y'r(a“v,5“-“V)-en(ii)dii . 

Without loss of generality assume |q;i| < A^/2. From Lemma 3, the last line is 

<Csup\ [ \d^^f{t,x,v)fdv] ' ||5“-“VWII- 
a:GR3 Id J 

Again, we first use the embedding 1F^’®(K.^) C L°°(R^) and second the L®(R^) 
Sobolev inequality to see that the above is 

0<|a|<2 0<|a|<2 

The last inequality follows from (10). This completes the estimate of h{f). □ 

Next we estimate the electric field E{t,x) in terms of f{t,x,v) through the 
macroscopic equation (19) using Lemma 5. 

Lemma 6. Say [f,E,B] is a classical solution to (4) and (7) with (8). Let the 
small amplitude assumption (10) he valid for some Mq < 1. 3(7 > 0 such that 

y E WEit)\\<\\ii-p)f\\+ E \\9vm- 

a|<Af-l 0<|Q|<Af 

Lemma 5 is used to remove the hydrodynamic part, P/, from this upper bound. 
Proof. We use the plus part of the macroscopic equation (19): 

-2d^Ei = + a“h+ - d^^dtb, - a“a*a+. 

By (15), 

Wd'^dM + ||a“9*a+|| < ciWPd^^dtfW + ||P9“97II} < c{||5“9j|| + ||a“97ll}- 

Since Mq < 1 in assumption (10), applying Lemma 4 to d'^l^- and Lemma 5 to 
we deduce Lemma 6. □ 
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We now prove the positivity of (5) for a classical solution of small amplitude to 
the full Vlasov-Maxwell-Boltzmann system, [f{t,x,v),E{t,x),B{t,x)]. The proof 
makes use of the estimates established in this section (Lemma 4, 5, 6) in addition 
to an exact cancellation property of the macroscopic equations. 


Proof of Theorem 2. From Lemma 1 we have 

It thus suffices to show that if (10) is valid for some small Mq > 0, then there are 
constants Ci,C 2 > 0 such that 

^ \\p^‘^fit)\\l<c^ m-p)dvmi + c 2 ^. 

0<|a|<Af |a|<Af 


Here G{t) is defined in Theorem 2. By (15), we clearly have 

i||P9“/(t)||2 < P“[a+,a_]|p + + ||a“c|p. 

The rest of the proof is therefore devoted to establishing 


^ {||a“K,a_]|p + lia“6|p + P“c|p} < c \\ii-p)dyit)\\l 

0<|a|<iV |a|<Af 

(26) +CMo 

0<\a\<N 

dG(t) 


+Co- 


dt 


which implies Theorem 2 when Mq is sufficiently small. This is because the last 
term on the right can be neglected for Mq sufficiently small: 

Wf{tW = ||P5“/(f)|p + ll(i-p)a“/(t)lP 

< {||a“[a+,a_]|| + ||a“5|| + P“c||r + 11(1 - P)9“/(t)|p. 


To prove (26), we estimate the macroscopic equations (16) through (20). 


The estimate for b{t, x) with at least one spatial derivative. We first esti¬ 
mate Vxd°‘b with |q;| < fV — 1. We take of (17) and (18) to get 




3 


Y d^^d°^hi I + 


Yi-d^^d^bj + &d^kj + + {d^d^k + 

3¥=i 


dtd^d'^c} 


3¥=i 

Pd^d°‘li P d^d°‘h^ - dtd^d°‘c. 

Since dtd^d^^c = 2i9ti9*9“c, we get 

= dtd^d'^c + - d^d°‘hj + d^d°^hj + d^d°^hij} + d^d°^h + d^d°^hi 

3i^i 


= + d^d^kj + d^)9“/iy} -b 2{d^d°‘U + 
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Therefore, multiplying with yields: 

(27) ||V,9“6.|| < cY,{\\d-l^\\ + \\d-h,\\ + \\d‘^k,\\ + \\^^ 

This is bounded by the first two terms on the r.h.s. of (26) by Lemma 4 and Lemma 
5. We treat the pure temporal derivatives of b{t,x) at the end the proof. 

The estimate for c{t,x). From (16) and (17), for |q;| < TV — 1, we have 

||V,a“c|| < P%|| + P“/le||, 

||5ta“c|| < p*5“&,|| + ||a“T,|| + ||a“/i,||. 

By (27), Lemma 4 and 5 both ||9t9“c|p and ||Va;9“c|p are bounded by the first 
two terms on the r.h.s. of (26). 


The estimate for [a+{t,x),a-{t,x)]. By (20), for |a| < N — 1, we get 
\\dtd‘^[a+,a.]\\<\\dVt\\ + \\d‘^h^\\. 

By Lemma 4 and 5, ||9t9“[a+, a_]|| is thus bounded by the first two terms on the 
r.h.s. of (26). 

We now turn to purely spatial derivatives of a±{t, x). Let |a| < — 1 and 

a = [0, oi, 02, as] ¥= 0- 

By taking 9* of (19) and summing over i, we get 

(28) -A,9“a± ± V, • d^E = V, • dtd^b - ^ + h^}. 

I 

From the Maxwell system in (21) we have 

V,,-9“F; = 9“a+-9“a_. 

We plug this into (28), then multiply (28) with d°‘a± respectively and integrate 
over Adding the two resulting equations yields 

||V,9“a+||^ + ||V,9“a_||2 + ||9“a+ - d^a.\\^ 

= f ^9*9“a± dx. 

i,± 

We therefore conclude that 


||V,9“a+|| + ||V,,9“a_|| < ||9*9“6|| + ^ + h^}\\. 

i,± 

Since a is purely spatial, this is bounded by the first two terms on the right side of 
(26) by (27), Lemma 4 and 5. 

Now we consider the case 0 = 0. The same procedure yields 

i||V.a+||2 + i||V,a_||2 


We redistribute the derivative in t, as in [12], to estimate the first term by 


/ d 
WX ■ dtb{a++ a-)dx = — 

d 

< — 

dt 


VX ■ b{a+ + a-)dx — J Vx ■ bdt{a+ + a_) 

V,, • b{a+ + a_)dx + C\ |V^fejp + C EllftOilP. 

± 
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The first term on the right is dG/dt and the last two terms have already been 
estimated by the other terms on the right side of (26). 

The estimate for b{t, x) with purely temporal derivatives. Lastly, we con¬ 
sider dtd°'bi{t, x) with a = [oq) 0, 0,0] and oq < N—1. In [11], Guo used a procedure 
which involved the Poincare inequality and Maxwell’s equations (21) to estimate 
d‘^dtbi from the macroscopic equation (19), which includes the electric field. But 
the Poincare inequality does not hold in the whole space case to estimate the low 
order pure temporal derivatives of b. 

On the other hand, there are two equations for dtbi in (19). We can add them 
together to cancel the most dangerous term, Ei, and obtain 

2dthi + d^a+ -{■ d^a- = I'^i + h'^i + l,,^ + 

Without this exact cancellation, our estimates for Ei are not sufficient to conclude 
the required estimates for dtd°^bi{t,x). From this equation we can estimate the 
temporal derivatives of b directly: 

\\did‘^h{t,xW < cY, (P*9“a±||2 + + h^}\\^) . 

± 

If do > 0 then the terms on the right have already been estimated above by the 
first two terms terms on the right side of (26). 

Alternatively if oq = 0, we use the previous estimate for d''a± to obtain 

\\dMt,x)\\^ <Cj^ J V, • 6(a+ + a_)dx + q|V,&||2 

+ C^(||5*a±|p + ||{Z± + 4}|H. 

± 

From Lemma 4 and 5, we finish Theorem 2 with G = Vx • b(a+ + a-)dx. □ 

5. Global Existence 

The main goal of this section is to establish global existence of classical solutions 
to the Vlasov-Maxwell-Boltzmann system, (4) and (7) with (8). Our proof uses 
Theorem 2 and energy estimates (Theorem 3) developed in this section. 

We have the following global existence theorem: 

Theorem 3. Let [f{t, x, v), E{t, x), B{t, x)] be a classical solution to to the Vlasov- 
Maxwell-Boltzmann system (4) and (7) with (8) satisfying (10). Then there exists 
an instant energy functional satisfying 

j^£{t)+v{t) < ./mm- 

With (10), this yields global existence from a standard continuity argument. 

First we prove the following nonlinear estimate: 

Lemma 7. Let jaj -I- |/3| < N. Then there is an instant energy functional such that 

I (a^r(/, /), 9|(i - p)/) I < £^^/t)v{t). 

The main improvement of this estimate over Lemma 3 is the presence of the 
dissipation, I?(t). Glearly ||i?|p is not needed in the dissipation in Lemma 7. 
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Proof. First assume (3 = 0. We have 

9“r(/,/)= ^ c“r(9“-“v,a“V)- 

cni <a 

Furthermore (see for instance [6, p.59-60]): 

(29) (r( 5 ,p),r) = (r( 5 ,p),(I-P)r). 

Plugging these last two observations into Lemma 3 yields 

|(a“r(/,/),a“/)| < c ^ (|a“-“VI.|5“Vl2 + |a“-“Vl 2 |a“VI.) I(i- p)5“/|.. 

ai <Q; 

Further integrate over to obtain 

|(a“r(/,/),9“(i-p)/)| 

<cV / (|a“-“vi.|5“vi2 + |a“-“vi2|a“vi.)l(i-p)5“/|.dx. 


Q;i <Q; 


By Cauchy-Schwartz and symmetry the above is 


<c||(i-p)a“/iu ^1/ (|5“-“vi.|a“Vl2)'dx 

ai<a 


1/2 


Our goal is to show that 


Jk3 J 

This will imply Lemma 7 in the case (3 = 0. We will use the embedding 
L°°(R®) followed by the L®(]R®) Sobolev inequality to one of the terms inside the 
X integration. But, because the dissipation does not contain Pf, it is not enough 
to just take the supremum in x over the term with the smaller number of total 
derivatives. Instead, if jo; — ai| < JV/2 then we always take the sup on the i/ norm: 

f (|5“-“VI.|5“Vl2)"dxj ' < ||5“VII sup 

1/2 


c 


< c'ii3“vii '^{v)\\d°‘ “v(v)iUi.6(R3)C^v 

<C'I|5“VII y" I / v{v)\d°‘-‘^^d^V^f\^dxdv \ ' < 

|a|<l ^•^«^XR3 J 

Alternatively, if \a — ai\ > N/2 then |q;i| < N/2 and f\\i, is part of the 

dissipation so that it is ok to take the supremum on the other norm to obtain 

I [ (|a“-“VI.|5“Vl2)'dxj ' < ||a“-“^/|U sup |5“V(a;)|2 
ljR3 J a;eR3 

< C'||5“-“VIU yy I / \d°‘'^d^V^f\^dxdv\ ' < 

UR3xR3 j 

This completes the proof of Lemma 7 when /3 = 0. 

If \(3\ > 0 then we expand / as in (11) to obtain 

r(/, /) = F(P/, Pf) + r(p/, (I - P)/) + r((i - p)/, p/) + f((i - p)/, (i - p)/). 
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We will estimate each of these terms using similar Sobolev embedding arguments to 
above. However, we use this splitting because 9/3?/ is not part of the dissipation. 
The exponential velocity decay of P/ is used to control the hard sphere weight (6). 
Alternatively (I — P)/ is part of the dissipation. We only estimate the first two 
terms below. The other two terms can be treated similarly. 

Lemma 3 and the definition (15) yield 

|( 9 |r(p/,p/), 9 |(i-p)/)| 

< C ^ [|9“-“^p/|2|9“ip/|2 + |9“-“^P/|2|9“^P/|2] |9|(I - P)/|.. 

ai <Q: 


Further integrate the above over use Cauchy-Schwartz and symmetry to obtain 

|( 9 |r(p/,p/), 9 ^(i-p)/)| 


<C||9|(I-P)/|U Y. \ [ 

ai<a 


1/2 


Without loss of generality assume |ai| < N/2. Then by the same embedding 



pf\l\d^^pf\ldx< 


sup |a“i/(a;)| 



P/IP 


\|a|<i / 

This yields the estimate in Lemma 7 for the dpT{Pf,Pf) term. 

Finally, we estimate the i9^r(P/, (I — P)/) term. In this case Lemma 3 yields 

|(9^r(p/, (I-P)/), 9^(1-p)/) I 
< C||9|(I - P)fhY {- P)/| 2 ) dxj 

+q|9|(I-P)/|U^|y^^(|9“-“^P/|2|9“ni-P)/l.) dx^ . 

Since (I — P)/ is always part of the dissipation, we can take the supremum in 
X on the term with the least number of total derivatives and use the embedding 
i7^(R.^) C to establish Lemma 7 for this term. And the last two terms can 

be estimated exactly the same as this last one. □ 


We will now introduce a bit more notation. For 0 < m < N, a reduced order 
instant energy functional satisfies 

Y \\d0fitw+ Y 

|/3|<m,|a| + |/3|<Ar |a|<Ar 

Similarly the reduced order dissipation rate is given by 

v^it)^\\Eit)r+ Y ii^“/wii'+ E ii9|{i-p}/(t)ip. 

0<|a|<A l/3|<m,|a| + |/3|<Ar 

Note that, £N{t) = £{t) and I?Ar(<) = 'D{t). With this notation in hand, we can 
now prove the main global existence theorem of this last section by induction. 
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Proof of Theorem 3. We will show that for any integer m with 0 < m < N we have 

(30) + V^{t) < 

Theorem 3 is a special case of (30) with m = N. We will prove (30) by an induction 
over m, the order of the v—derivatives. We may use different equivalent instant 
energy functionals on different lines without mention. 

For m = 0 and |a| < N, by taking the pure 9“ derivatives of (4) we obtain 

{dt + v-V, + qiE + vx B)- V49“/ - {d‘^E ■ vjyTJIq, + L{9“/} 

(31) = - ^ C^^q{d°‘^E + vxd°^^B)-Wvd°^-°‘^f 

ai^O 

+ H {|(9“if;- v)9“-“v + r(a“V9“-“v)}. 

cti <a 

We will take the T^(R^ x inner product of this with 9“/. The first term is 
clearly 5^||9“/(t)|p. From the Maxwell system (7) we have 

- {d^E ■ vy^qi,d^f) = -d^E ■ d^J = ~ {Wd^EitW + P“S(t)||2) . 

Next up, since v{v) > c(l + |u|), we use two applications of Cauchy-Schwartz to get 

[ \q{d°^^E + vx d°‘^B) ■ V„9“”“4| \d°‘f\dxdv 

< [ (|9“if;| + |9“iB|)| [ i/(v) |v„9“-“4pdul ' [ [ iyiv)\d°‘ff dv\ ' dx 

Jr3 (7r3 J UR3 J 

< ||9“/|U (|9“^F;| + |9“^i?|)" 44 I V.9“-“V|' ' . 


If |a| = 0 this term is non-existent. Therefore, for this term, 0 < |a| < and 
|q;i| >0 so that either |q;i| < iV/2 or |a — aij -|-1 < iV/2. Without loss of generality 
say |ai| < iV/2; we use the embedding i7^(R^) C L°°(K^) to obtain 

sup (|9 “if;| -k |9“iB|) < C V (||9“9“iS(t)|| -k ||9“9“iB(t)||) < a/^. 

|a |<2 

We can put the last two estimates together and use |a| > 0 to obtain 


/R3xR3 


Iqid’^^E + vx d°‘^B) ■ V„9“-“V |9“/| dxdv 


< ii9“/ii.v^iiv„9“-“4iu < ^/mm- 


It remains to estimate the last two nonlinear terms on the right in (31). 

First we estimate |(9“iif • v)d°‘~°‘^ f. Consider a = 0 which implies oi = 0. We 
bound for ||if|| by the dissipation to estimate this case. We split / as in (11) so 
that the exponential velocity decay of the hydrodynamic part controls the extra |u| 
growth: 


{f{r-v)u) 


< c 


/R3xR3 


\E\\v\{\Pf\^ + \{l-P)f\^)dxdv. 
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By the Sobolev embedding C we have 

[ \E{t,x)\\v\\iI-P)ffdxdv<c(snp |£;(t,a:)|) ||(I-P)/||^ 

jR3xR3 \£ceR3 / 

^ l|a“s(t)|| I \\{1-P)f\\l<^)v{t). 

Vl«l<2 / 

For the term with Pf, we obtain 

[ \E\\v\\Pf\^dxdv < C f sup \Pfix)\ 2 ) [ \E\\Pf\ 2 dx. 
JR3xR 3 \a:GR3 / JR3 

By the Cauchy-Schwartz inequality 

f \E\\Pf\ 2 dx < q|F;(t)||||p/|| < 


We just controlled ||Fl(t)|| by the dissipation, For the other term, from 

(15), the embedding C L°°(]R^) and the Sobolev inequality we obtain 


(32) 


sup |P/(x)| 

,a;GR3 


< C||[a,6,c]|Ui,e(K3) <CY, llV.a“/|| < CI?i/2(t). 

|a|<l 


Adding together the last few estimates we obtain 


(|{E ..)/■/) 


< 


C [ \E\\v\{\Pf\^ + \{l-P)f\^)dxdv<^/6{i)V{t). 

jR3xR3 


This completes the estimate when a = 0. 

Alternatively if |a| > 0, we use iy{v) > c(l + |w|) and Cauchy-Schwartz twice as 

<C [ A| p(v) \d°‘f\ dvdx 

^R3xR3 

<C f \d‘^^E\\ [ dv\ ' I [ iy{v)\d°‘ffdv\ ' dx 

JR3 (Jr3 J (Jr3 J 

< c\\d‘^f\\:. ff dv^ dx^'. 

Here ||9“/||,y is controlled by the dissipation since |a| > 0. But we could also have 
|a| = |ai| = -/V and then would not be controlled by the dissipation. Therefore 
this estimate requires more care. 

If |a — ai| < lV/2, then to ensure that we have at least one derivative on / we 
use the embedding IF^’®(]R^) C followed by Sobolev’s inequality to obtain 


(l{d^^E-v)d°^-°^^f, a“/j 


sup [ p(u)|a“ °‘^f{x,v)\‘^dv<cf p(T)||a“ “V(u)||^i,6(R3)dw 

£!;eR3 JR3 JR3 

< C ^ ||V,9“-“19“/||2 < CV{t). 

|a|<l 

Alternatively, if |q; — ai| > N/2 then |ai| < N/2 and 

(33) sup \d‘^^E{t,x)f < C||9“^F;||^i,a(R3) < C V < C6{t). 

a:eR3 
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In either case the last two estimates establish 


\d°‘^E\‘ 


IR3 


v{v) 9“ “V dv>dx< S{t)'D{t). 


Combining these last few estimates, since |q;| > 0, we obtain 

< s/WW)- 

It remains to estimate the nonlinear collision operator in (31). 

But from Lemma 7 and the basic invariant property (29) we have 


Adding up all the estimates for the terms in (31) and summing over |q;| < iV yields 


(34) 


l_d 
2 dt 


E l|5“/f+ ||9“[A,i3]f 

|q|<A 


+ E (i9“/(t),9“/(t)) < ^)V{t). 

\a\<N 


By Theorem 2 and Lemma 6 for some small dg > 0 we see that 


2^1 E Wff + W\E,B]f-2C^G{t)\+5’Mt)<^)V{t). 
[|a|<Ar J 

We add (34) multiplied by a large constant, Cq > 0, to the last line to obtain 
(35) j^Soit) + W) < 

We have used L > 0 from Lemma 1. Above, for Cq large enough and from the 
definition of G{t) in Theorem 2, the following is a reduced instant energy functional: 


Soit) 



We conclude (30) for /3 = 0. 

Say Theorem 2 holds for |/3| = m > 0. For |/3| = m + 1, take 9^ of (4) to obtain: 

{dt+vV^ + q{E + vx B)- V49^(I - P)/ - 9“£; • dpiv^jqi 

+9^{L9“(I - P)/} + E 
01^0 

(36) = E |{9“^ A • 9;3,u}9“I“;/ + 9^r(/, /) 

- E E C“C|^g9;3,ux9“^il-V.9“r“;/ 

lail^^O |ai| + |/3lM0 

-{dt + v-V^+qiE + vx B)- V49^P/. 


Notice that we have split the equation in terms of microscopic and hydrodynamic 
parts in a different way from (12). We do this to estimate the linear part of the 
collision operator properly in terms of our weak whole space dissipation. Standard 
estimates for (9/3^9“/, 9^/), e.g. Lemma 2, are too strong in this case. 
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We take the inner product of (36) with 9^(1 — P)/ over x We estimate 
each term from left to right. The first inner product on the left gives 

i|lia|(i-p)/(t)|p. 

Since \(3\ = m + 1 >0, |a| < N — m — 1. Then by Lemma 6 the second inner 
product on the left is bounded by 

-1 {d^E ■ d0{v^}q^,d$(l - P)/) I > -CWE\\ ||5|(I - P)/|j 

>-c^\\d^Ef-T^\\d^pii-v)fr 

>-C,Po(t)-?7l|5|(I-P)/f. 

For the linear operator, since \dp^v(v)\ is bounded for /3i ^ 0 (Lemma 2) we have 
(a^{p5“(i - p)/}, 9^(1 - p)/) = Wp(i - p)f\\l 
+ ^ cl (I - P)/, 9|(I - P)/) 

> \\d^{l-P)f\\l-C\\d^(l-P)f\\ ||5|_^^(I-P)/||. 

/3i5^0 

Then together with Lemma 2, since L = v — K, we deduce that for any rj > 0 there 
is a constant > 0 such that the third term on the left side of (36) is bounded 
from below as 

(5;3{L9“(I-P)/},a^(I-P)/) > ||5|(I-P)/||2-,7 ^ \\d^,il-P)f\\l 

\0'M0\ 

-c, Y. P|(i-p)/ip. 
m<\0\ 

We can further choose Crj > 0 such that the inner product of the last term on left 
side of (36) is bounded by 

77lia|(i-p)/(t)|p + c, ^ 

l/3i|=l 

We split the last term above as in (11) to obtain 

<r;P^(I-P)/(t)||2 + C, ^ \\W,d^_^Pf\\^+C^ Y l|V.a^_^,(I-P)/|P 

l/3i|=l l/3i|=l 

<7j\\d^(i-p)fm‘^+cvo{t)+c, Y iiv.5|_^,(i-p)/ir 

l/3i|=l 

This splitting was used to get rid of the velocity derivatives of the hydrodynamic 
part, which are not in the dissipation. 

Next we estimate all the terms on the right side of (36). Let us first consider 
the first term on the right, using (1 + |i;|) < ciy{v) we have 

9|(i - P)/) 

<c [ ,yiv) dT$if \d^{i- p)/| dxdv. 

jR3xR3 


(37) 



20 


ROBERT M. STRAIN 


We use Cauchy-Schwartz twice to see that the last line is 

<C [ [ v{v) ' \ [ ,.{v)\d^{l-P)f\^dv] ' dx 

Jr3 IJk3 ^ ^ J IJr3 J 

We split / into it’s microscopic and hydrodynamic parts (11) so that we can remove 
the velocity derivatives from the hydrodynamic part: 




+qia|(i-p)/lj.|y^j9“^i?r|y^^K^’)|a;i;,ni-p)/| dv'^dx'^ . 

Since the total number of derivatives is at most N, either |q;i| < N/2 or |q; — q;i|-|- 
\(d — (di\ < N/2. So we take the supremum in x of the term has the least total 
derivatives, for the microscopic term, using the embedding C to 

observe 

f \d^^E\^U i,iv)\d;z^^{i-p)ffdv\dx<sit)v{t). 

J ]R^ L J 

On the other hand, for the hydrodynamic variables we have 

f \d'^^E\^\f i^iv) dv\dx<C f E^ \d°^-'^^Pf\ldx. 

Jr3 [Jr3 ^ ^ J Jr3 

If |a — Oil < N/2, then as in (32) we see that 

[ |a“-“ip/|2dx < C\\d‘^^E\\'^V{t) < C£{t)V{t). 

Jr3 

If alternatively |a — q;i| > N/2 and |q;i| < N/2 then we use (33) to get 

[ \d°^^Ef |a“-“ip/|2dx < C'£:(t)||5“-“VII^ < C£it)V{t). 

Jr3 


In any case, we conclude 

- p)f) I < Vmm- 

This is the estimate for the first term on the right side of (36). 

For the nonlinear collision operator in (36), by Lemma 7 we have 

|(a^r(/,/),9|(i-p)/)| < Vmm- 

For the next two terms on the right side of (36) we have 

(qd^^E ■ V„9“-“V + qd 0 ,v X 9^(1 - P)/) 

<C [ v{v)(\d<^^E\ V„9“-“V +\d^^B\ |a^(I-P)/|dudx. 

Jr3xR 3 \ ^ H H y 

Notice that these are in the form (37). They are therefore estimated the same way. 
Furthermore the same procedure (as in the estimate for (37)) yields 

I {qiE + vxB)- V45|P/, 9|(I - P)/) I < Vmm- 
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It remains to estimate 

\{{dt + vv,}d^pf,d^ii-p)f)\. 

Since |/3| > 0, we have |a| < TV — 1 and the above is 

<C [ iy{v) {\dtd^Pf\ + |V,9|P/|) |9|(I - P)/| dxdv 

dR3xR3 

< c\\d^(i - p)/||. {Hd^Pfll + l!v,9|P/|U) 

< c\\d^{i - p)f\\A\\dtd‘^Pf\\ + l|v,9“P/||) 

<7j\\d^ii-p)f\\l + c^ ^ Wd^^dvr- 

\a\ = l 

This completes all the estimates for the terms in (36). 

By collecting all the estimates for the terms in (36) and summing over \f3\ = m+1 
and |a| + |/3| < iV we obtain 

E { - cMt) 

|/3|=m+l,|a| + |/3|<JV ^ 

<Cv E \\d^{l-P)fit)\\l 

|/3|=m+l,|a| + |/3|<Af 

+c, E \\d^ii-p)f{m" + Vmm- 

\d\<m,\a\ + \/3\<N 

Here C is a large constant which does not depend on rj. Choosing rj = ^ we have 

E U jtW - p)/wf + w - p)/wii^| 

|/31=m+l,|a| + |/3|<Af ^ 

(38) < CVm{t) + yW)T^it)- 

Now add the inequality from (30) for |/3| = m multiplied by a suitably large con¬ 
stant, Cm, to (38) to obtain 


jACmSm{t) + l E \\d^iI-P)f{t)r\ +{Cm-C)Vm{t) 

[ |/3|=Tn+L|a|-S|/3|<JV J 

E \\d^{i-p)fmi<^/W)m- 

\f3\=m+l,\a\ + \/3\<N 

Here Cm is chosen so that Cm — C > 0. We further define 

Sm+lit)=Cm£mit) + ^ E || (I - P)/(0 f • 

|/5|—m+l,|a| + |/3|<A^ 

Since ||9/3P/|j < C|jP/||, £m+iit) is an instant energy functional. □ 
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